New decay modes of gravitino dark matter 
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We consider the three-body decays of gravitino dark matter in supersymmetric scenarios with 
bilinear i?-parity violation. In particular, gravitino decays into £W* (iff') and vZ* (yff) are 
examined for gravitino masses below Mw- After computing the gravitino decay rates into these 
three-body final states and studying their dependence on supersymmetric parameters, we find that 
these new decay modes are often more important than the two-body decay, into a photon and a 
neutrino, considered in previous works. Consequently, the gravitino lifetime and its branching ratios 
are substantially modified, with important implications for the indirect detection of gravitino dark 
matter. 



INTRODUCTION 

Recent measurements indicate that about 25% of the 
energy density of the Universe is made up of a mysterious 
form of non-baryonic matter known as cold dark matter 
[lj]. To explain it, physics beyond the Standard Model 
is required. Supersymmetric extensions of the Standard 
Model are the most promising scenarios that may account 
for the dark matter. They are motivated, in addition, by 
the hierarchy problem and by the unification of the gauge 
couplings. 

Within supersymmetric models, R-parity conservation 
was long believed to be a prerequisite for supersymmetric 
dark matter, for it guarantees the stability of the lightest 
supersymmetric particle (LSP) -the would-be dark mat- 
ter candidate. It was pointed out in Q, however, that 
if the gravitino -the superpartner of the graviton that 
arises in local supersymmetric theories- is the LSP, it 
can be a suitable dark matter candidate in R-parity vio- 
lating models. Indeed, even though the gravitino is un- 
stable in such scenarios, its lifetime may be much longer 
than the age of the Universe. Such a long lifetime is the 
result of the gravitino feeble interactions, which are sup- 
pressed by the Planck scale and by small R-parity violat- 
ing couplings. A scenario with gravitino dark matter and 
R-parity violation is even favoured by thermal leptogen- 
esis, as it alleviates the tension between the large reheat- 
ing temperatures required by leptogenesis and the con- 
straints from Big-Bang Nucleosynthesis In R-parity 
violating models, therefore, the primordial gravitino pro- 
duced in the early Universe is a viable and well-motivated 
dark matter candidate. 

An important feature of these R-parity violating mod- 
els is that, unlike their R-parity conserving counterparts, 
gravitino dark matter can be indirectly detected. In fact, 
the small fraction of gravitinos that have decayed until 
today constitutes a source of high energy cosmic rays 
d, In their decays,_gravitinos may produce gamma 
rays [1, Q , neutrinos 0] and antimatter Q that could 



be observed in present and future experiments. Present 
data from FERMI, for instance, already constrains the 
parameter space of these models [9] . 

The indirect detection signatures of gravitino dark 
matter strongly depend on the gravitino lifetime and on 
its branching ratios. In previous works, it has been as- 
sumed that they are determined by the two-body decays 
of the gravitino. We will show, in this paper, that that is 
not always the case. For gravitino masses below Myy, the 
three-body final states £W* and i/Z*, where * denotes a 
virtual particle, typically give a large contribution to the 
gravitino decay width, modifying in a considerable way 
the gravitino lifetime and its branching ratios. These 
results have important implications for scenarios of de- 
caying gravitino dark matter, especially regarding their 
indirect detection prospects. 



TWO-BODY DECAYS 

We work in the framework of supersymmetric models 
with bilinear R-parity violation with gravitino as LSP, 
such as those considered in [1, H, @]. In these models, 
the sneutrino field acquires a vacuum expectation value 
that breaks lepton number and gives rise to mixing be- 
tween neutralinos and neutrinos and between charginos 
and charged leptons. Gravitinos may then decay into 
standard model particles via lepton number violating in- 
teractions. 

Because the R-parity violating couplings are expected 
to be largest for the third generation, we will assume, 
following [7j , that the sneutrino acquires a vev only along 
the i> T direction. It is convenient to define a dimensionless 
R-parity breaking parameter as £ r = (y T )jv, where v = 
174 GeV. This parameter is constrained from neutrino 
masses and from successful BBN to be in the range [5j 
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In this setup, all R-parity violating processes, including 
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the decay of the gravitino, are controlled by £ r . 

Gravitinos can decay into two-body final states con- 
taining either a neutrino (and a photon, Z boson, or 
higgs boson) or a charged lepton (and a W boson). Be- 
cause R-parity is broken along the v r direction, the final 
state charged lepton is a r and the neutrino is a v T . The 
gravitino decay rates can be calculated directly from the 
interaction Lagrangian [l(| • For the two-body final states 
7^ T , Zv T and W + t~ , they are given by [H, 0] 
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and P, f,j, h are simple functions -see Q for their def- 
inition. In this expression, are ' res P ec ~ 
tively, the photino-zino, the zino-zino, and the wino-wino 
mixing parameters Ml. To a good approximation they 
can be written as frl Till] 
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where sw and cw denote the sine and cosine of the 
weak mixing angle, and M\,M 2 are the U(l) and SU(2) 
gaugino masses. For simplicity, we will assume that 
gaugino masses are universal at the GUT scale, so that 
M 2 ~ 1.9Mi at the electroweak scale. Notice, from the 
above equations, that the mixing parameters decrease 
with gaugino masses. 

The two-body decays of gravitino dark matter have 
been studied in detail in previous works 0, EH . For 
gravitino masses below Mw, gravitinos decay into ^v T 
with a 100% branching ratio. Above the W threshold, 
this decay mode becomes quickly negligible, accounting 
for less than 1% of the branching for gravitino masses 
larger than 150 GeV. If heavier than the W boson, the 
gravitino dominantly decays into W ± t^ , with significant 
contributions from Zv T (once it is open) and, in some 
cases, from hv T . The main point of this paper is that to 
properly compute the gravitino lifetime and its branch- 
ing ratios the two-body decays considered so far in the 
literature are not enough; a certain class of three-body 
decays must also be taken into account. 
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FIG. 1: The feynman diagrams contributing to the gravitino 
decay into rW* . 



THREE-BODY DECAYS 

For gravitino masses smaller than Mw the only two- 
body channel kinematically available for the decay of the 
gravitino is 7^ T . In addition to it, there exists three-body 
final states consisting of a lepton and a virtual massive 
gauge boson, such as tW* (t//') and v T Z* {v T ff), that 
may contribute significantly to the gravitino decay width. 
These three-body final states are considered for the first 
time in this paper. 

Two diagrams contribute to the decay of the gravitino 
into the three-body final state tW* (rff) - see figure 
[T] The contribution of the (a) diagram is proportional 
to Uww ' wn ereas (b) , which comes from the non-abelian 
4-vertex, is independent of it. We have computed, from 
these two diagrams, the decay rate T(G — > t~ + W + * — > 
t~ //'). The differential decay rate can be written as 

rfr 



dsdt 256tt 3 7 
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where N c is the color factor, and s and t are, respectively, 
the invariant masses for the //' and the fr systems. The 
squared amplitude can be factorized as 
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where / is a long function of the given parameters. As 
usual, the decay rate is obtained from the differential one 
by integrating over the possible values of s and t. 

An important feature of the function / is that, due 
to the interference between the diagrams (a) and (b) in 
figure [1] it contains terms quadratic in, linear in, and in- 
dependent of Uyy^y. These latter terms, which are anal- 
ogous to those found in equation ([J}, will play a crucial 
role in our analysis, as they favour the decay into the 
three-body final states for larger gaugino masses. 

The gravitino decay into v T Z* proceeds through dia- 
grams analogous to those from figure [1] and the expres- 
sion for its decay rate has a similar form. In principle, 
the decay into v T ^* should also be taken into account, as 
it interferes with the v r Z* contribution. We have veri- 
fied, however, that the photon-mediated diagram is very 
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FIG. 2: The ratio between the gravitino lifetime computed 
from two- and three-body decays and that one obtained from 
the two-body final state only as a function of the gravitino 
mass for different values of Mi. 



much suppressed and can be neglected. In bilinear R- 
parity breaking models, other three-body final states are 
not expected to be relevant for gravitino decays. 

The total decay width of a gravitino with m,Q < M\y 
will, then, be given by 

Ftot 2 (G) - r(G -> 7 ^)+r(G -»• t+w-*)+y{g -> v T z*), 

(ii) 

where the factor two takes into account the charge con- 
jugated final states. This decay width is proportional to 
the R-parity breaking parameter see equation (|10|). 
The gravitino lifetime, rg, is simply given by the inverse 
of the decay width, l/T tot (G). In the following we study 
the gravitino lifetime and branching ratios as a function 
of supersymmetric parameters. 

To begin with, let us consider the effect of the three- 
body final states, tW* and v T Z* , on the gravitino life- 
time. Figure [5] displays, as a function of the gravitino 
mass, the ratio between the correct gravitino lifetime 
(computed from two- and three-body final states) and 
the gravitino lifetime obtained from the two-body final 
state only. This ratio, which is independent of £ T , would 
be equal to 1 if the contribution from three-body final 
states were negligible. Clearly, that is not the case. The 
lines in that figure correspond to different values of the 
gaugino masses. As expected, the heavier the gauginos, 
the larger the effect of the three-body final state. No- 
tice that there is a significant deviation from the two 
body result over a wide range of gravitino and gaugino 
masses. For Mi — 300 GeV, for instance, the effect of 
the three-body final states is significant already for grav- 
itino masses around 40 GeV, yielding a lifetime about a 
factor two smaller than that computed for two-body final 
states; as the gravitino mass increases the effect becomes 
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FIG. 3: The gravitino branching ratios as a function of the 
gravitino mass for Mi = 100 GeV. 



larger, reaching an order of magnitude close to Mw- For 
Mi = 1 TeV, the genuine gravitino lifetime can be more 
than two orders of magnitude smaller than that obtained 
from two-body decays. Undoubtedly, the three-body fi- 
nal states tW* and vZ* have a considerable impact on 
the computation of the gravitino lifetime. 

Besides the lifetime, the other important quantities 
that determine the indirect detection signatures of grav- 
itino dark matter are its decay branching ratios. Figure[3] 
compares the two-body and three-body decay branching 
ratios of the gravitino as a function of m,Q for Mi = 100 
GeV. For this value of Mi the dominant decay mode of 
the gravitino is the two-body final state "fv T , but with a 
branching ratio that can be much smaller than 1, reach- 
ing about 60% for mg ~ Mw- We also notice from the 
figure that, among the two three-body decay final states, 
it is the tW* that dominates. Its branching can reach 
30% for ?7ig ~ Mw- Remarkably, even for small gaugino 
masses, the contribution of the three-body final states is 
significant. 

For slightly larger gaugino masses, the branching ratios 
look quite different, as illustrated by figure |H This figure 
is analogous to the previous one but for Mi — 300 GeV. 
Notice that, in this case, the two body final state domi- 
nates only for gravitino masses below 50 GeV. From that 
point up to Mw, it is the tW* final state that accounts 
for most gravitino decays. Close to Mw, the tW* final 
state has a branching of about 70%. It is also observed 
in the figure that the vZ* branching can be quite signif- 
icant, reaching 20% and becoming even more important 
that the two-body decay for gravitino masses between 70 
GeV and Mw- This figure demonstrates that three-body 
final states may easily dominate the gravitino decays. 

To better study its dependence with gaugino masses, 
we show, in figure O the branching ratio into three-body 
final states, BR(G -)• tW*) + BR(G -> v T Z*), as a 
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FIG. 4: The gravitino branching ratios as a function of the 
gravitino mass for Mi = 300 GeV. 



function of the gravitino mass for different values of M\ . 
Notice that the three-body final states are always signif- 
icant and that they dominate the gravitino decays even 
for moderate values of the gaugino masses. If M\ = 1 
TeV, the three-body final states have a branching larger 
than 90% over the entire gravitino mass range we con- 
sider, 40GeV < mg < M w . 

In view of these new results on the gravitino lifetime 
and its branching ratios, existing constraints on scenar- 
ios with decaying gravitino dark matter have to be reex- 
amined, and the prospects for the indirect detection of 
gravitino dark matter have to revised. Because of the 
smaller branching into 7^ TJ the gamma ray line will be 
less pronounced than previously believed, and it will be 
accompanied by a continuum of gamma rays coming from 
the three-body decays of the gravitino. Besides, the pos- 
sibility of discovering gravitino dark matter via antimat- 
ter opens up for gravitino masses below M\y, thanks to 
the positrons and antiprotons generated by the tW* and 
i> T Z* final states. In a future work, we will study these 
interesting topics in detail. 

CONCLUSIONS 

We have studied the decays into tW* and v T Z* of 
gravitino dark matter in the context of supersymmetric 
models with bilinear R-parity violation. We computed 
the gravitino decay rates and showed that these previ- 
ously neglected processes typically give large contribu- 
tions to the gravitino decay width. They actually dom- 
inate the decay rate over a wide range of gravitino and 
gaugino masses. In a future work, we will study in de- 
tail the implications of these results for gravitino dark 
matter. 
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FIG. 5: The gravitino branching ratio into three-body final 
states as a function of the gravitino mass for different values 
of Mi. 
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